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Key points 

• Basic models are developed for active lava flow surface disruption at a slope break 

• Surface disruption can significantly affect core cooling and flow length, depth and 
advance rate 

• Disruption length scales depend on flow regime (laminar vs turbulent) beyond slope 
break 
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Abstract Changes in the underlying slope of a lava flow impart a significant fraction of 
rotational energy beyond the slope break. The eddies, circulation and vortices caused by this 
rotational energy can disrupt the flow surface, having a significant impact on heat loss and thus 
the distance the flow can travel. A basic mechanics model is used to compute the rotational 
energy caused by a slope change. The gain in rotational energy is deposited into an eddy of 
radius R whose energy is dissipated as it travels downstream. A model of eddy friction with the 
ambient lava is used to compute the time-rate of energy dissipation. The key parameter of the 
dissipation rate is shown to be pi? 2 ///, where p is the lava density and u is the viscosity, which 
can vary by orders of magnitude for different flows. The potential spatial disruption of the lava 
flow surface is investigated by introducing steady-state models for the main flow beyond the 
steepening slope break. One model applies to slow-moving flows with both gravity and pressure 
as the driving forces. The other model applies to fast-moving, low-viscosity, turbulent flows. 
These models provide the flow velocity that establishes the downstream transport distance of 
disrupting eddies before they dissipate. The potential influence of slope breaks is discussed in 
connection with field studies of lava flows from the 1801 Hualalai and 1823 Keaiwa Kilauea, 
Hawaii, and 2004 Etna eruptions. 

Index terms 8425 Effusive volcanism 

8414 Eruption mechanisms and flow emplacement 
8485 Remote sensing of volcanoes 
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1. Introduction 


When a flow changes direction, the momentum vector of the flow also changes direction. 
This results in the production of angular momentum when the flow transits a break in slope or is 
directed laterally by confining topography. Under relatively steady eruption conditions, a viscous 
lava flow rapidly forms a solid crust that insulates the interior of the flow. If undisturbed, this 
crust preserves the mobility of the hot lava core and allows the flow to travel substantial 
distances until the lava supply ceases or the interior of the flow cools enough to inhibit further 
advance. A comprehensive review of the influence of cooling, crust fonnation and effusion rate 
is contained in Harris and Rowland [2009]. 

The angular momentum generated by slope breaks and lateral confinements results in some 
form of circulation, eddying or vortex formation within the flow or at its margins. Pulses in lava 
supply can also induce circulation in the flow. When such flow patterns disturb the cooler or 
crusted surface of the flow, they can have a profound influence on radiative heat loss, viscosity, 
incipient crystallization and thus the distance a flow travels [ Finch and Macdonald, 1953; Booth 
and Self, 1973; Moore, 1987; Rowland and Walker, 1990; Crisp and Baloga, 1990a, b; Crisp and 
Baloga, 1994; Harris and Rowland, 2001]. 

The only predictive attempt at calculating the extent of surface disruption appears in Harris 
and Rowland [2001]. It is based on empirical studies of core exposure as a function of flow 
velocities in channels and they pointed out that this dependence lacked rigorous quantification 
[Harris and Rowland, 2001, p.28]. Field measurements and theoretical modeling of the main 
flow from the well-documented 1984 Mauna Loa eruption show a marked increase in viscosity 5 
km downstream of a major slope break at 10 km along the flow path [Moore, 1987; Harris and 
Rowland [2001]. With the slope embedded in the computation of flow velocity, the FLOWGO 
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model of Harris and Rowland [2001] predicts such an increase [see Figure 5 in Harris and 
Rowland, 2001] and motivates further study of slope changes and surface disruption. 

There have now been more than two decades of remote sensing observations of active lava 
flows in numerous settings. Such images commonly show relatively broad spikes in spectral 
radiance (and derived metrics such as radiant flux) along the flow path that are indicative of 
transient or persistent disruptions of the cooler flow surface and a consequent exposure of the 
hotter inner core [James et al., 2007; James et al., 2010]. Wright and Flynn [2003] and Wright 
et al. [2010] give methods for estimating the core exposure of active flows from high resolution 
satellite remote sensing data. Although such data show indications of potential relationships 
between slope changes and the disruption of the lava surface, higher spatial resolution remote 
sensing data are needed to definitively identify a cause and effect relationship, and to place 
observational constraints on the length scales of disruption. 

The work presented here attempts to predict the extent of flow surface disruption and core 
exposure due to a sudden change in the underlying slope. The model is a combination of three 
first-order models for different aspects of the problem. A basic mechanics model is used to 
approximate the rotational energy caused by a steepening slope change. Similar energy 
considerations apply to a shallowing slope change, lateral redirection of the flow by topographic 
confinements, channel constrictions and wall effects, breaches, and pulses in the lava supply rate; 
the resultant impact of each case on the surface disruption may be different. The gain in 
rotational energy is deposited into an eddy whose energy dissipates as it is transported 
downstream by the main flow. An elementary model of the friction between the eddy and the 
ambient lava provides the rate of energy dissipation as a function of time. The downstream 
spatial disruption of the flow surface is subsequently determined with steady-state models for the 
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main flow beyond the slope break: one for slow-moving flows driven by both gravity and 
pressure, and another for low-viscosity turbulent flows. Both models of the main flow include a 
volume conservation requirement at all points along the path of the flow. 

The influence of surface disruption on cooling of the lava core is investigated. Implications 
in the context of surface disruption are given for field studies of basaltic flows from eruptions at 
Mt. Etna, Sicily in 2004, and the Hualalai 1801 and 1823 Kilauea eruptions in Hawaii. 

2. Model Overview 

The intent of the model is to develop a first-order estimate for the distance over which the 
surface of a steady-state flow is disrupted by a steepening of the underlying slope from one 
inclined plane to another. When a surface disruption occurs, it persists for a finite time as it is 
transported downstream by the underlying main flow. This duration is referred to here as the 
‘disruption time’. The corresponding ‘disruption length’ is obtained from the disruption time by 
knowing the flow velocity. The disruption length is determined by how much rotational energy is 
imparted by the slope break and the assumption that this energy goes into vertical circulation 
within the flow. This rotational energy then dissipates with time due to friction with the ambient 
lava. The model assumes that disruption time is given by the dissipation time. 

The model considers a vertical column of height h, width w, and length dx along the flow 
direction. The width and density of the flow are taken as constants. The column encounters the 
break in slope, then rotates to a new orientation at the beginning of the downstream segment 
(Figure 1). The pivot of the column is taken as the slope break and no slippage is considered. 
For illustration purposes, the slope of the upper reach is taken as zero and the lower reach has 
slope <f> with respect to the horizontal. The assumption is made that the rotational kinetic energy 
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imparted by the slope break is determined by the amount of reorientation of the column to the 
nonnal of the new slope. 

As a flow moves downslope, there is of course a change in kinetic and potential energy. In 
steady state on a constant slope, the kinetic energy of a column remains constant, but the energy 
gained from the decrease in potential energy is identically dissipated by the viscous stresses and 
is converted to heat. The disruption model presented here investigates the additional rotational 
kinetic and potential energies above the steady state main flow values. In absolute terms the 
energy changes at the slope break could be very small compared to the downstream steady state 
main flow values. However, the transfer of material from the upper layers of the inner core 
resulting from the additional rotational energy can have a significant influence on the thennal 
balance. 

The conservation of kinetic, potential and rotational energy is used to find the new boundary 
conditions at the beginning of the downstream segment of the flow. The volumetric flow rate is 
conserved throughout the upstream and downstream segments of the flow. Two cases are 
considered for the nature of the flow on the lower slope. One considers the flow to be essentially 
laminar with the possibility that both gravity and pressure are the driving forces. The other uses 
a gravity-driven hydraulic model more appropriate to a turbulent lava flow. Various details of the 
main flow models appear in the appendices. 

The induced circulations, eddies, and vortices are modeled approximately as a rotating 
cylinder oriented transverse to the main flow direction. Although many types of flow 
circulations are possible, they are referred to collectively as ‘eddies’ throughout this work. The 
model assumes that the eddies cause a surface disruption that persists essentially as long as do 
the eddies. The eddies dissipate within the flow as they are transported downstream. In the 
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model, the eddies generated by the slope break are simply carried by the main flow and draw no 
energy from the change in downstream potential energy. The dissipation rate is computed by 
estimating the energy lost due to eddy circulation against a viscous ambient fluid. This lost 
energy is removed from the eddy as the main flow travels downstream until the energy loss 
accumulates to a final value. For a given eddy radius, angular velocity, and viscosity, the 
persistence time can thus be calculated and determines the downstream eddy propagation 
distance. The variables for the quantitative description of this model are given under Notation. 

2.1 Mechanics at the Slope Break 

In going from one slope to the other, the normal to the original upstream surface rotates 
through an angle 0(t) from 0 at t = 0 to a new position with (Xtj) = </>. At the completion of the 
rotation, the column is aligned with the new normal to the downstream surface (see Figure 1). 
The time it takes to complete this rotation is tf. At t = 0, just before rotation begins as the column 
encounters the slope break, the column has translational kinetic energy, 

KE = Vi mu 2 , (1) 

where m is the mass of the column and u is the flow velocity upstream of the slope break. 

Across the slope break the column is treated as a physical pendulum of length h and angular 
velocity <x>. The center of mass of the column drops due to the rotation about the pivot point 
adding potential energy to the column in the amount 

PE = Vi mgh( 1 -cosr/9 , (2) 
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where g is the acceleration due to gravity. The sum of the kinetic energy and the change in 
potential energy goes into two quantities, rotational kinetic energy and translational kinetic 
energy immediately after the slope break. 

Using the approximation for the moment of inertia of the column, / = mh 2 /3, the angular 
acceleration is given by 

dd_ 3gsmO(t) ' wherg 0 (O ) = 0, co = co(0) = u/h . (3) 

dr 2h 

The elementary solution where sin# ~ 6 is 


6(t) = 




< < p , where 



(4) 


The initial rotational kinetic energy, KE rot , at the slope break is given by 


1 1 mh 2 

KE rot (0) = ~Ico(0f =- — 


(uX 


V 


1 2 

= —mu , 
6 


h) 

for h > 0 . 


(5) 


Equation (5) shows that as the column begins to pivot over the slope break, the rotational 
energy is one-third that of the incident translational kinetic energy given by equation (1). At the 
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completion of the rotation about the pivot point, at t = tj, the rotational kinetic energy imparted 
by the slope break is 


KE rot {t f ) = -Ico{t f Y = 


2 1 mh ( U | / tslT -tjr 


2 3 2 h 


e f + e J 


= ~ mu 2 1 e ‘ fh + e ‘ fh j , for h > 0 . 


( 6 ) 


Across the slope break an energy balance can thus be constructed. The incoming kinetic 
energy from equation (1), plus the change in rotational kinetic energy from equations (5) and (6), 
plus the gain in potential energy from equation (2), must go into the outgoing kinetic energy of 
the column, 


mu 2 lco(0) 2 

~2T + 2 



Iu> 2 f 

~Y 



( 7 ) 


where (Of is the angular velocity of the column after it has rotated through 0 to its new 
orientation. Equation (7) with flow rate conservation ( uh = U2/12) can be solved to obtain the 
outgoing velocity, U 2 , on the immediate downstream side of the slope break with the selection of 
the physically appropriate roots. The same type of energy balance applies when a slope shallows. 
The main difference would be the subsequent boundary conditions and the nature of the 
governing transport equation beyond the slope break. A more refined approach would 
investigate a continuous deformation of the column across the slope break, but that requires 
solution of the Navier-Stokes equation and is beyond the scope of the present work. 
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The exit value of the velocity uj changes very little (<0.5%) from the incident conditions for 
slopes (<15°) relevant to most lava flows. This is because the energy balance in equation (7) 
predominantly puts the gains in potential energy into the rotational energy that increases with 
slope, leaving the exit flow conditions essentially unchanged. For practical purposes, therefore, 
these velocity changes can be ignored. 

The absolute value of the rotational energy per kg changes substantially with the slope and 
the incident conditions, as shown in Figure 2. The vertical axis gives the exit rotational energy 
KE 10 ,{tf) per kg as a function of the change in underlying slope. The relative changes in the 
rotational kinetic energy are illustrated by considering three flow depths (h = 1,3 and 6 m; gray, 
red, and green curves, respectively) and two incident flow velocities ( u = 1 and 2 m s , solid and 
dashed curves, respectively). The rotational energy acquired from the slope break increases with 
incident flow thickness and changes by factors of 2, 4, and 7 (gray, red, green curves) for the 
thicknesses shown over slope changes up to 15°. When the incident velocity is doubled (dashed 
curves), there is the same dependence on flow depth and slope, but there is an additional overall 
increase in the rotational energy, as expected. A much greater sensitivity to slope is found by 
changing the flow depth keeping the velocity constant. 

2.2 Eddy Dissipation 

Eddies are produced by the angular momentum generated by the slope break. In the first- 
order model used here, the eddy is treated as a cylinder of fluid with radius R that rotates as it is 
transported downstream by the main flow. The best case for the persistence of an eddy is when 
its surface is in contact only with the ambient fluid that transports it. In the worst case, an eddy 
contacts counter-rotating neighboring eddies and breaks into smaller ones. 
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At time t = tf, the rotation of the lava column due to the slope break is completed and an eddy 
is established. The time of eddy dissipation is then measured with time beginning at tf The 
rotational kinetic energy from the slope break is deposited into rotational energy of the eddy, 
thus 


KE {t r) = —I co 2 = — 

rot v j 7 ^ e eo ^ 



( 8 ) 


where m e and I e are the mass and moment of inertia of the eddy, respectively. The initial angular 
velocity of the eddy rotation, co eo , can be found by re-arranging equation (8) to be 


4KE rot (* f ) 2 

KE ro M f ) 2 

m KEJf f ) 

R 2 m e R ]] 

m e R 

m m 


( 9 ) 


where the rotational kinetic energy per unit mass is obtained from the energy balance for the 
column in equation (6). 

The resistive force on the eddy due to ambient lava is approximated by 


F-ftA — 
e dr 


t>A e co e (t ) , 


( 10 ) 


where A e is the surface area of the eddy in contact with ambient lava. Assuming the eddy retains 
its integrity, cumulative energy lost due to friction over time beginning with t/is given by 


W dm(t)=( Fds « uA (j f co e (t)ds = u(2jiRw)( co e (t)—dt 

J'f J'.r J, f at 


— at. t>t 


f 


( 11 ) 
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The variable s is the distance a point on the cylinder travels to dissipate the energy W,/, ss by 
friction with the ambient fluid, given by 



= Rc °eo~ U-e 


- 4 (/ 


-y )/r ) . 


( 12 ) 


To determine the amount of dissipation after the onset of eddy rotation, <x> e (t) must also be 
detennined. The time dependence for the deceleration of eddy rotation due to friction is found 
from Newton’s 2 nd law 


j 

e dt 


= -FR - -juA e co e R . 


(13) 


Integrating equation (13), we obtain 


co e (t) 

f 

OJ eo 


d(o e 

Me 



(14) 


and an expression for a> e {t) is found to be, 


co (t) = a) e 

e v ' eo 


-4(t-t f )/T ' 


where T = 


P# 


(15) 
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This shows that the critical parameter for eddy dissipation is the time constant T. It can also 
be expressed in terms of the Reynolds number of the flow itself, Re = puih^/p as T = R eR 2 /q 0 . 
Fast-moving, high-Re flows will have long time constants and the effects of slope breaks will 
tend to propagate greater distances downstream and vice versa. 

The rotational energy lost to friction is given by equation (11), using equations (12) and (15), 
as 


II 


_ g-w-i/Vr )_ KE m (i/) (l - e ' m ' ,,VY ), 


(16) 


which yields the time decay of rotational energy due to the friction of the eddy. Curves relevant 
to lava flows are shown in Figure 3. Slightly more than half (55%) the initial rotational energy is 
dissipated by friction in 0.1 time constant, 98% is dissipated in half a time constant, and 99.99% 
is dissipated when t-tf = T. Thus T is a reasonable estimate of the dissipation time and, by 
assumption, the disruption time. 

2.3.1 Eddy Size 

Many factors could influence the circulation of a flow as it traverses a change in slope or 
encounters a topographic barrier. These include the character of the incident flow (laminar, 
turbulent, or disrupted), roughness of the flow bed, the presence of entrained and incipient solids, 
and interactions with channel walls. All such factors require independent theoretical or empirical 
studies. 

There is a link between the volume of the fluid column and the volume of the eddy that is the 
recipient of the rotational energy imparted by the slope break. To embrace a range of eddy sizes, 
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the extent of the column along the direction of flow ( dx ) is prescribed to be the diameter of the 
eddy (2 R; Figure 1). Thus the fraction of column volume that receives the rotational energy is 
given by 

o, _ vol. eddy 
vol. col. 

For a large circulation, R might approach h/2 as suggested by engineering experience with 
turbulent flow in pipes. It is difficult to estimate a minimum relevant size of an eddy in a lava 
flow. However, from a practical standpoint a lower limit can be estimated by considering an 
incident laminar flow. Appendix A gives the cumulative rotational energy for an incident laminar 
flow as it traverses the slope break. Due to the strong dependence of rotational energy on the 
height within the flow, only 3% of the rotational energy would come from the lower half of the 
flow and only about 11% from the lower two-thirds. The upper 5% of the flow contributes 
almost a quarter of the rotational energy and so R = 0.1/7 is taken as a practical minimum value 
for eddy radius. 

The time constants for eddy dissipation shown in Table 1 are clearly dominated by the 
viscosity and the eddy size, which must be roughly proportional to the flow depth. Figure 3 
shows the relative dissipation as a function of time for T = 100, 200, 500 and 1000 s from 
equation (16). Low-viscosity lavas imply large time constants, and such flows are predicted to 
maintain circulations for long times, tens of seconds or more. Eddies would be carried substantial 
distances downstream by any fast-moving segment of the flow. If the flow is turbulent, however, 
the simplistic model used here would require significant modification to account for the decay of 
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larger eddies into smaller ones. Conversely, high-viscosity flows could manifest the influence of 
circulation by local cooling, stagnation, and inflation. 

The dominant role of viscosity in determining the extent of surface disruption calls into 
question the validity of viscosity estimates derived from field measurement of active flows. 
Usually such estimates are derived from the Jeffreys’ equation for a steady-state flow of constant 
depth [Nichols, 1939], although numerous modifications have been made to account for time- 
dependent effects [e.g., Baloga and Pieri, 1986], lava yield strength [e.g., Harris and Rowland, 
2001], levee building [e.g., Glaze et a/., 2009], and similar factors. These formulations do not 
account for the fluid pressure caused by the resulting topographic gradients in the flow. It is 
shown in Appendix B that the incorporation of seemingly small pressure gradients causes an 
important systematic departure from the steady-state constant depth assumptions of the Jeffreys’ 
approach. The primary consequence is that evaluations of viscosity based on field measurements 
of flow parameters could be overestimated by an order of magnitude or more in the case of thick 
(50-200 m) flows on other planetary surfaces. 

3. Implications 

3.1 Influence of Flow Surface Disruption on Core Temperature 

The potential influence of surface disruption on the core temperature is illustrated by 
considering a schematic flow of constant depth /?2 and velocity ic that begins with temperature T 0 
at the top of an inclined plane. The core temperature is computed by the basic steady state 
radiation loss fonnula 
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LoJx) = T 0 


. 3 f eoT] h,x 

1 H — — 2 — 

pc p U 2 


(18) 


using the flow parameters given in Table 2, with / denoting the areal fraction of exposed hot lava 
core within the crusted flow surface, and a being the Stefan-Boltzmann constant. 

The degree to which the disrupted upper layer of a lava flow exposes core lava depends upon 
many factors, including the viscosity of the lava and the thickness and mechanical strength of the 
crust. We do not attempt to draw a direct link between the rotational eddies and the ability to 
actually “disrupt” the surface crust. However, several cases exploring a range of potential / 
values can be explored to place bounds on the sensitivity of increased surface disruption over 
short distances. 

Figure 4 shows the resulting temperature profiles along the length of the flow, for three cases 
having different values of f The uppennost curve assumes a core exposure fraction of /= 0.05, 
which is somewhat high for most flows [ Crisp and Baloga, 1990a; Oppenheimer, 1991; Wright 
et ah, 2000; Wright et al., 2010]. For comparison, the lowennost curve uses an extreme value of 
/= 0.5. The red curve shows the strong influence of the 1 4 radiation term when the flow surface 
is disrupted by eddies. Twelve 60 m segments of disruption with /= 0.9 were inserted in the first 
2 km of the flow, followed by three others further downstream. Everywhere else, a value of / = 
0.05 was applied. The disruption in the first 2 km of the flow causes a major drop (35°C) in the 
core temperature. Using equation (3) of Harris and Rowland [2001], such a temperature drop 
would cause the core viscosity to increase by a factor of 4. This would have a significant 
influence on the advance velocity and the depth of the flow. A few segments of surface 
disruption upstream will have a significant impact on the ultimate flow length. 
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3.2 The September 2004 Lava Flow at Mt. Etna 

A relatively well-documented flow-producing eruption exemplifying laminar flow (Appendix 
B) occurred at Mt. Etna, Sicily in 2004-2005 [e.g., Mazzarini et al., 2005; Burton et al., 2005]. 
Comprehensive dimensional data on the flow were acquired by airborne laser altimetry 
Mazzarini et al. [2005]. Analysis of 162 profiles indicated a typical flow thic kn ess of 5 m, 
although there is considerable variability, with a maximum of 17 m. Both Mazzarini et al. [2005] 
and Wright et al. [2010] show that 20-30° slopes are common, and field observations obtained 
from the Global Volcanism Network [2004] suggest that 1 m s' 1 is a representative flow velocity. 
The eddy dissipation relationship (Figure 3) shows that the influence of slope changes (or 
equivalent influences) decays in -0.5-2 min. for T ~ 100-200 s. With channel velocities of 1 m 
s' 1 , this corresponds to disruption lengths of 30-120 m. In Appendix B, it is shown that the use of 
the Jeffreys’ equation (commonly used for flow velocity calculations) can significantly 
overestimate the viscosity of active flows, so a dissipation time on the order of hundreds of 
seconds may be reasonable. If the lava viscosity was significantly in excess of hundreds of Pa s, 
the dissipation time could be on the order of seconds to a few tens of seconds for eddies of < 1 m 
and a flow thickness of 5 m. However, significant slope changes occur about every 100 m along 
the 2004 Mt. Etna flow path. Even with a relatively rapid dissipation rate, the repeated slope 
changes at length scales comparable to the dissipation lengths keep adding flow circulation and 
surface disruption. 

One potentially useful approach to validating and placing constraints on the proposed model 
is to correlate changes in lava flow surface temperature, derived from remote sensing thermal 
imaging, with underlying surface topography and slopes. However, the spatial resolution of even 
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some of the best orbiting thermal imagers, such as Hyperion (30 m pixels), is insufficient. In 
order to distinguish changes in surface radiance on the order of the disruption scales (30 - 120 m 
for a flow like the 2004 Mt. Etna flow) requires spatial resolutions about an order of magnitude 
finer (~5-10 m pixels). One approach for future work is to use field thermal imagers that are 
capable of spatial resolutions on the order of 5-10 m/pixel (or better), with specific application to 
this problem to ensure synoptic coverage over the full length of possible disruption. Such field 
campaigns have successfully documented variations in lava flux [James et al., 2007; James et 
a/., 2010]. In particular, James et al. [2007] examine the time dependent thermal flux of several 
lava flows during the 2004-2005 eruption at Mt. Etna that progress from a steeper to a shallower 
slope. A similar technique would be appropriate here for lava flowing onto a steeper slope. It is 
important to note that such a field campaign must adequately characterize the time dependent 
changes in thermal flux owing to pulses in lava flux in order to distinguish the systematic surface 
disruption due to the change in slope. 

3.3 The 1801 Hualalai Flow, Hawaii 

The possibility of turbulent lava flows has been postulated for many years [e.g., Nichols, 
1939, p.294; Shaw and Swanson, 1970; McGetchin and Eichelberger, 1975; Baloga et al., 1995]. 
Turbulence can be caused not only by high flow velocities and low viscosities, but also 
constrictions and widenings of the flow, small and large-scale topographic variations, and the 
motions of entrained crystals and ambient materials. Although there is some disagreement on the 
eruption rates during the 1801 Hualalai lava flow [Kauahikaua et al., 2002], with viscosity 
determined by petrologic studies to lie in the range 1 0 1 0 2 Pa s, and field estimates of flow 
velocities of ~10 m s' 1 , the 1801 Hualalai eruption likely produced at least transient turbulent 
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lava flows with depths of ~5 m near the source. The uncommonly deep (6-18 m) downstream 
channels may have formed in part through construction due to very short duration overflow 
events, however, even partially full, greater turbulence through these stretches is expected 
[McGetchin and Eichelberger, 1975; Guest et al, 1995]. 

The 1801 Hualalai flow also traversed several large changes in slope (3-6°) along the flow 
path [Baloga et al., 1995]. Xenoliths were deposited as bedload at the beginning of the final 
reach toward the ocean. From that point, the decay of these circulations and intense cooling 
produced a morphology common to other basaltic channelized flows on low slopes. 

The elementary turbulent flow model in Appendix B indicates that such a flow would reach 
the terminal velocity on a flat plane in tens to -100 m from the source or slope break. Using the 
range of likely terminal velocities (5-15 m s' 1 ) suggested in Baloga et al. [1995], the computed 
disruption lengths are shown in Figure 5 for appropriate Hualalai parameters (p = 2600 kg m' 3 , h 
= 5 m, R\ = 1 m, (j) = 10°) and viscosities of 50-500 Pa s spanning the range identified by 
McGetchin and Eichelberger [1975]. With the eddy radius taken as only 20% of the flow depth, 
the computed flow disruptions extend for hundreds of meters to a kilometer or so. The 
circulations gained from the slope breaks most likely contributed significantly to the suspension 
of the xenoliths. Given that there are several slope breaks along the flow, such disruption lengths 
are consistent with the complex morphologic features of the flow [ McGetchin and Eichelberger, 
1975; Baloga et al., 1995; Guest et al., 1995] until the very shallow slopes near the ocean. 

3.4 The 1823 Keaiwa “Great Crack” Flow, Kilauea Volcano, Hawaii 

The 1823 Keaiwa flow from the Great Crack fissure at Kilauea also produced rapidly moving 
low -viscosity flows. Near the source Guest et al. [1995] estimated the velocity at 15 m s' 1 . Field 


19 


420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 


evidence suggests a source depth of <0.5 m with a downstream thickening in distal regions to 
only 1-2 m near the ocean [Guest et a/., 1995; Ba/oga et al., 1995; and sources cited therein]. 
The predominant underlying slope is about 5° from the fissure to the ocean (a distance of -4 km 
using Area 2 of Figure 4 from Baloga et al. [1995]). The flow apparently issued from the Great 
Crack as a rapidly moving sheet, transitioning to a slabby ‘a’ a at the distal margins. The 
turbulent main flow model (Appendix B), using the parameters u( 0) = 15 m s' 1 , h( 0) = 0.5 m, (fh 
= 5° and h-i = 2 m, indicates that the flow would attain a velocity U 2 of 3.75 m s' 1 in -100 m. 
Thus, the 4 km transit time to the ocean would be -1000 s. The model used in Baloga et al. 
[1995] assumed a constant velocity of 10 m s' 1 from the fissure to the ocean resulting in a transit 
time about half this value. With a velocity of 3.75 m s' 1 , R - 0.5-1 m and // = 100 Pa s, Table 1 
suggests a disruption length of about 24-98 m. Due to the thinness of the flow, small-scale 
topographic variations could occur over such a disruption length, feeding the slope induced 
circulations. The resulting mixing would have contributed to the homogeneity of the advancing 
sheet, but the transit time was evidently too short to permit a significant increase in the viscosity. 

4. Summary and Conclusions 

Changes in slope or the lateral redirection of a lava flow impart a significant fraction of the 
incident kinetic energy of the flow into rotational energy. For steepening slopes, the eddies, 
circulation and vortices caused by this rotational energy can disrupt the flow surface and have a 
significant impact on the heat loss and thus the distance the flow can travel. The quantity of 
rotational energy imparted to the downstream flow is more sensitive to the flow depth than the 
incident flow velocity. There is a relatively large quantity of potential energy available from the 
main flow compared to the rotational energy of eddies that disrupt the surface and alter the 
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thermal balance. The preliminary model presented in this work does not draw energy from the 
main flow. Repeated slope changes and rough topography could significantly extend the 
disruption of the flow surface. Many lava flows experience multiple slope breaks of various 
magnitudes that can combine to disrupt the surface. A more refined theoretical treatment and 
comprehensive field measurements are needed to explore this possibility. 

This work provides only a first-order analysis of the essential steps in disrupting a lava flow 
surface to the extent that it could affect the thermal heat balance of the flow. Consequently, there 
are numerous opportunities for improvement at each step of the analysis. First, in transitioning an 
abrupt slope break, a continuum approach based on Navier-Stokes equations would provide a 
refinement. Relaxing the assumption of an abrupt slope change is also a more realistic approach. 

The physics of eddy dissipation is another area for future study. Circulations, particularly if 
turbulence is generated, are three-dimensional, unlike the planar assumptions used here. This 
suggests that a horizontal component of rotation and dissipation must also be included. Lateral 
confinements, wall effects, and changes in the lateral direction of the flow path could contribute 
significantly to the rotational energy of the flow and thus to surface disruption. Unlike the simple 
estimate used here, the physical sizes of the eddies diminishes as they propagate downstream. 
The details of how these considerations affect the heat balance are reserved for a future study. 
Nonetheless, such improvements in the theory might provide a basis for future field studies and 
approaches to flow and eruption dynamics inferred by remote sensing. 

Besides the parameters of the incident flow, this work shows that the primary factors 
controlling the surface disruption are the size of the eddies and the viscosity of the downstream 
ambient lava. The model results obtained in this work are at least qualitatively consistent with 
surface disruptions interpreted for the fast-moving flows from the 1801 Hualalai and 1823 
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Keaiwa eruptions, as well as a recent flow at Mt. Etna. Future avenues for developing model 
constraints include comparison of predicted disruption lengths with high spatial resolution 
thennal remote sensing data. Future modeling work will include investigation of other likely 
contributors to circulations and surface disruption in active lava flows. 


Appendix A: Cumulative Rotational Kinetic Energy 

Here the cumulative rotational kinetic energy of the column in laminar flow as it encounters an 
abrupt slope break is calculated. Within the flow, the velocity profile is taken as 


u{z) = 


pg sin 6 { 

p 



(Al) 


where hi is an arbitrary height within the flow interior up to the top surface at h. The upstream 
surface is assumed to be inclined to the horizontal at an angle 0\. If there is no inclination of the 
upper surface, the flow is driven by a constant pressure (momentum flux) and the calculation 
procedes replacing gsin#i with gdh/dx 
The rotational kinetic energy up to hi is 


KE (/?,) = 

rot v 1 ' 


I(h x )a>{h x ) 2 1 phpvdx 

2 2 3 


1 ph\wdx 

gsin^/zj 2 

2 3 

1 2 ^, j 
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where, v is the lava kinematic viscosity (= p/p). Note that the moment of inertia and the angular 
velocity depend on hi and the strong dependence on the flow depth. Up to the full height h of the 
flow we have: 


KEJh h ) = 


1 pg\sm6 { ) 2 h 5 wdx 

24 V 1 


(A3) 


Thus the cumulative rotational KE fraction as a function o f h \ is 

r.mmjhd. ( A4) 

KE(h) [ h ) 

This shows, as expected, that the very upper layers of a laminar flow contribute the vast 
majority of rotational kinetic energy as the flow goes over an abrupt slope break (see Figure Al). 

Appendix B: Main Flow Transport 

The main flow modeling referred to in the text is presented below. The analysis of the 
mechanics of the slope break provides the boundary conditions for main flow transport of the 
eddies beyond the slope break. The depth and velocity of the main flow may change as a 
function of distance beyond the slope break, but the volumetric flow rate is conserved at all 
locations downstream. 

The laminar case 
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The governing equation for the flow depth is taken as 


dh 2 g sin (phi dh 2 

dt v dx 


gcostp d 
3v dx 


/ 3 dM 

/?2 5 

V dx ) 


(Bl) 


where v is the kinematic viscosity of the lava (= /-i/p). The second term on the left-hand side 
represents gravity as a driving force while the right-hand side represents the influence of fluid 
pressure. The pressure term is often ignored in volcanologic applications, but is a critical term in 
hydraulic formulations. Here we consider the steady-state solution of equation (Bl). The first 
integration gives 


q 0 =u 2 { 0)h 2 (0) = 


g sin (phi 
3v 


g cos (p 3 dh 2 N 
3v 2 dx i 


u 2 (x)h 2 (x) , 


(B2) 


where q Q is the volumetric flow rate per unit flow width. A more convenient fonn is given by, 


1 = 


(h,\ 

\ h J> 


1 - cot (p 2 


dh 2 

dx 


hj = ? I 


3 7A 
g sin (p. 


(B3) 


where hj is the oft-cited Jeffreys’ steady-state flow depth. 

Three types of possible behaviors are evident from equation (B3), depending on the boundary 
conditions (^(O), /z 2 (0)), the slope, and the viscosity. The different modes of flow behavior are 
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associated with three roots of equation (B2) or (B3). The critical value that determines the 
behavior of the flow is the Jeffreys’ flow depth hj. 

What is remarkable about the inclusion of pressure in the steady state is that the Jeffreys’ 
steady-state flow depth can only be attained for one fortuitous set of conditions of flow rate, 
slope, and viscosity and boundary conditions (7/2(0), h 2 ( 0)). For h 2 < hj (u 2 > uj,) the flow must 
thin and accelerate as it moves from x = 0. In steady state, the flow regime actually moves 
further away from Jeffreys’ conditions with distance from the source until the steady-state 
assumption is no longer valid, the viscosity increases due to cooling, or a traveling wave is 
established [Mei, 1966]. 

For h 2 >hj , the flow must thicken and slow. In this case, the pressure term will always 
attempt to drive the flow toward a pond-like topography with dh 2 ldx ~ sinf/n. This has nothing to 
do with whether there is a slope break upstream, it is solely a function of the driving forces. The 
regimes associated with equation (B3) are clearly analogous to supercritical and subcritical 
hydraulics concepts. The tendency to produce an almost horizontal upper surface topography on 
shallow slopes may be one of the principal factors causing pahoehoe lobe inflation [Hon et al., 
1994; Keszthelyi et al., 1999; Glaze and Baloga, 2013]. 

The analytic solution of equation (B3) is found by changing variables and the boundary 
condition 


ip (2f(x)) = h 2 (x) / hj , g(x) = tan (p 2 x / h r 
ip(0) = h 2 (0)l hj, £(0) = 0 


(B4) 


With these changes, the solution of equation (B3) requires integration of 
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r ip 3 dip _ 
J ip 3 - l 


(B5) 


The key to the integration of equation (B5) is the expansion of the denominator in terms of the 
real and two conjugate roots and the subsequent use of partial fraction expansions. 
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h- 1 ) 


i iS)( i Wfi 


-l 


lp + — + 

2 2 


1 i' 

ip + 

2 2 


ty 3 dip , 


\ 


(B6) 


where z'=V(- 1). With the use of the partial fraction expansions and considerable algebra, we find 


I 


ip 3 dz 

ip 3 - 1 


= ip + ln(ip - 1) 1/3 + tan 


V3 J y{3 


\ 


1 + 2 ip 


- In (tp 2 +ip + 1) 
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(B7) 


Thus the solution of equation (B3) in the ip variable is 


% = ip -ip{ 0) + ln 


( V'-l ] 

1/3 r 

+ ^ tan" 1 

f V3 ^ 

- 75 tan- 1 

f ^ 1 

-In 

' ip 2 +ip + l 3 
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(v'(O) 2 +^( 0) +1 j 


The analytic solution is found by undoing the change of variables from ip back to h using 
equation (B4). The analytic solution then gives the longitudinal flow profile h{x). Because the 
flow rate is constant along the flow path, equation (B2) then provides u(x). 
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Figure B 1 illustrates the dramatic influence of pressure on the longitudinal thickness profiles 
compared to a constant Jeffreys solution, for a flow rate and viscosity and slope that give hj= 5 
m. For clarity of illustration, the slope is taken as 0.2°. Once the flow rate, slope and viscosity 
are fixed, small increases in the flow depth at the source cause a significant departure from the 
constant Jeffreys’ flow depth. The most important consequence of this analysis is that field 
estimates using Jeffreys’ equation can significantly overestimate the viscosity. Figure B1 shows 
that at 3-4 km from the source, the flow has more than doubled beyond h j. Because viscosity 
estimates go as h 3 , the viscosity estimate would be overestimated by at least an order of 
magnitude. This sensitivity suggests that field measurements should measure the topographic 
gradient as well as the flow depth as is commonly done in hydrologic applications. 

The turbulent case 

Under conditions of turbulent flow, following the formalism of basic hydraulics, the re- 
oriented column on the lower reach (see Figure 1) experiences only two forces, gravity and flow 
bed resistance characterized by a single parameter, the bed stress. Consistent with a first-order 
analysis, the governing equation is taken as 

m C — L = m gs'm ( ft_o h A , (B8) 

dt 

where Ob is the stress at the flow bed and A is the area of bed contact with the flow. 

The common hydraulic application of this fonnulation is to empirically relate the bed stress 
under many different types of ambient conditions (e.g., bed roughness, channel geometry) to an 
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equilibrium flow depth. Similarly, a prescribed value of hi* is used here. With such an 
assumption, and the conservation of flow rate, equation (B8) reduces to 


du . . , 

u 2 = g sm< P 

ax 
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7 * \ 

fljr 
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with q = /z (0)u (0) . 
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The velocity is given by 


u(0)-u In 
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1 - au 2 
l-au(0) 2 


ag sin (px, 


with a = h* / q . 

2 J o 


(B10) 


This solution in general also has two regimes of flow behavior, but they approach the 
equilibrium flow depth hi* from above and below this depth. As the flow depth changes with 
distance, there are corresponding changes in the flow velocity to keep the flow rate constant. The 
simple model above assumes the flow is planar (i.e., infinitely wide), unlike actual flows with 
lateral confinement. Recent advances in hydraulic engineering show how to extend this 
elementary formulation to account for secondary flows within the main flow and the influence of 
narrow and wide stress-inducing lateral boundaries [ Guo andJulien, 2005]. 


Notation 

A area of contact between turbulent flow and flow bed (constant), nr 

A e area of contact between eddy and ambient lava (constant), m 2 

C p specific heat capacity of lava, J kg' 1 K' 1 

dx length of column in direction of incident flow (variable), m 
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/ fraction of exposed lava core 

F resistive force acting on eddy, N 

g acceleration due to gravity, m s' 2 

h depth of the incident flow (constant), height of column entering slope break, m 

hj equilibrium turbulent flow depth on lower surface (constant), m 

/? 2 flow depth on lower (2 nd ) surface (space dependent), m 

hi(0) flow depth on lower (2 nd ) surface after column rotation (constant), m 

hj Jeffreys’ steady-state flow depth, m 

/ moment of inertia of column (constant), kg m 2 

I e moment of inertia of eddy (constant), kg m 2 

KE translational kinetic energy of the incident flow (constant), J 

KE ro ,( 0) initial rotational kinetic energy of the incident flow (constant), J 

m mass of the column (constant), kg 

m e mass of eddy (constant), kg 

PE potential energy gained by column rotation (variable), J 

q a lava volume flow rate per unit channel width, m 2 s' 1 

r cumulative rotational kinetic energy fraction for laminar flow 

R radius of circulating eddy (constant), m 

s distance on the eddy through which the resistive force F acts (time dependent), m 

t time, s 

tf time required to rotate column through an angle (j) (variable), s 

T 0 initial lava core temperature (K) 

T core temperature of lava core (K) 

u velocity of the incident flow (constant), m s' 1 

U 2 flow velocity on lower (2 nd ) surface (space dependent), m s' 1 

i/ 2 ( 0 ) flow velocity on lower (2 nd ) surface after column rotation (constant), m s' 1 

w width of flow (constant), m 
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Wdi SS energy lost from eddy due to friction with ambient lava (time dependent) 
r time constant for eddy dissipation (constant), s 

s lava emissivity 

6 angle through which column rotates (time dependent) 

H dynamic viscosity of lava (constant), Pa s 

v kinematic, viscosity of lava (constant), m 2 s' 1 

p lava density 

o Stefan-Boltzmann constant, 5.670373(2 1 )x 1 0 s W m 2 K -4 

C p specific heat of lava 

T lava core temperature 

Ob bed stress for resistance to turbulent flow (constant) 

r time constant for angular rotation of column (constant), s 

<j) slope of the flow bed after the slope break (constant) 

co angular velocity of column (time-dependent) 

co e angular velocity of eddy rotation (time-dependent) 

co eo initial angular velocity of eddy rotation (constant) 

(Of angular velocity of column after it has rotated through angle <f> (constant) 

2s ratio of eddy volume to column volume 
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Table 1 . Dissipation time constants T in seconds for eddy of radius R and viscosity p. 


Eddy radius, 

Lava dynamic viscosity, p (Pa s) 

R( m) 

10 

100 

1000 

10000 

0.5 

65 

6.5 

0.65 

0.065 

1 

260 

26 

2.6 

0.26 

2 

1040 

104 

10.4 

1.04 

5 

6500 

650 

65 

6.5 

10 

26000 

2600 

260 

26 


Table 2. Parameters used to calculate lava core cooling 


Flow depth, h 
Initial temperature, T a 
Lava density, p 
Lava specific heat, C p 
Lava emissivity, e 
Flow velocity, u 


5 m 
1330 K 
2600 kg m 3 
1225 J kg' 1 K' 1 
1 

0.2 m s’ 1 
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Figure 1 . Geometry of a lava flow encountering a slope break. A lava column of height h and 
thickness dx (dashed lines) rotates through an angle 6 with angular velocity co as it pivots over 
the slope break, thereby imparting rotational energy to the flow and producing an eddy of radius 
R and angular velocity co e . On the new slope <p, the flow takes on new velocity u 2 and depth h 2 . 
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Figure 2. The absolute value of the exit rotational energy KE rot (tj ) per kg as a function of the 
change in underlying slope at the slope break. The reference case is taken as u\ = 1 m s' 1 and h\ 
= 1 m for the incident flow. The rotational energy changes by a factor of 2 over slope changes up 
to 15° (solid gray curve). For flows of thickness 3 and 6 m, the rotational energy changes by 
factors of 4 and 7 (red and green curves), respectively, demonstrating the sensitivity of rotational 
energy to slope for thicker flows. For comparison, the dashed curves show the effect of a greater 
incident velocity ( u\ = 2 m s' 1 ): the dependence on slope is essentially the same as for the lower 
velocity, but there is an overall increase in the rotational energy. 
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Figure 3. Eddy energy dissipation as a function of time for T = 100, 200, 500 and 1000 s from 
equation (16). 
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Figure 4. Theoretical temperature profiles along the length of a hypothetical flow for three 
different values of / The dark curves assume constant core exposure fractions of/= 0.05 and / = 
0.5. For comparison, the red curve shows the strong influence of the 7 4 radiation term when the 
flow surface is disrupted by eddies. Twelve 60 m segments of disruption with / = 0.9 were 
inserted in the first 2 km of the flow, followed by three others further downstream. 
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785 Figure 5. Disruption lengths as a function of viscosity for the lava flow from the 1801 Hualalai 

786 eruption for three plausible flow velocities. See text for discussion of parameters. 
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790 Cumulative rotational kinetic energy, KEJLhJ/KEJlh) 

791 Figure Al. Cumulative rotational kinetic energy for a laminar flow pivoting over an abrupt slope 

792 break as a function of the relative depth h\ within a flow of thickness h. 

793 
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Figure Bl. Comparison of steady-state gravity and pressure driven solutions with constant 
Jeffreys’ solution hj for the flow depth. For clarity of illustration, the slope is taken as 0.2°. 
Once the flow rate, slope and viscosity are fixed, small increases in the flow depth at the source 
cause a significant departure from the constant Jeffreys’ flow depth. At 3-4 km from the source, 
the flow has more than doubled beyond hj. 
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